Search for: All records

Gene regulatory divergence is thought to play an important role in adaptation, yet its extent and underlying mechanisms remain largely elusive for local adaptation with gene flow. Local adaptation is widespread in marine species despite generally high connectivity and is often associated with tightly linked genomic architectures, such as chromosomal inversions. To investigate gene regulatory evolution under gene flow and the role of inversions associated with local adaptation to a steep thermal gradient, we generated RNA-seq data from Atlantic silversides (Menidia menidia) from two locally adapted populations and their F1 hybrids, reared under two temperatures. We found substantial divergence in gene expression and thermal plasticity between populations, with up to 31% of genes being differentially expressed. Reduced thermal plasticity, temperature-dependent gene misexpression, and the disruption of coexpression networks in hybrids point toward a role of regulatory incompatibilities in local adaptation, particularly under colder temperatures. Chromosomal inversions show an accumulation of regulatory incompatibilities but are not consistently enriched for differentially expressed genes. Together, these results suggest that gene regulation can diverge substantially among populations despite gene flow, partly due to the accumulation of temperature-dependent regulatory incompatibilities within inversions.

 

Abstract The levels and distribution of standing genetic variation in a genome can provide a wealth of insights about the adaptive potential, demographic history, and genome structure of a population or species. As structural variants are increasingly associated with traits important for adaptation and speciation, investigating both sequence and structural variation is essential for wholly tapping this potential. Using a combination of shotgun sequencing, 10x Genomics linked reads and proximity-ligation data (Chicago and Hi-C), we produced and annotated a chromosome-level genome assembly for the Atlantic silverside (Menidia menidia)—an established ecological model for studying the phenotypic effects of natural and artificial selection—and examined patterns of genomic variation across two individuals sampled from different populations with divergent local adaptations. Levels of diversity varied substantially across each chromosome, consistently being highly elevated near the ends (presumably near telomeric regions) and dipping to near zero around putative centromeres. Overall, our estimate of the genome-wide average heterozygosity in the Atlantic silverside is among the highest reported for a fish, or any vertebrate (1.32–1.76% depending on inference method and sample). Furthermore, we also found extreme levels of structural variation, affecting ∼23% of the total genome sequence, including multiple large inversions (> 1 Mb and up to 12.6 Mb) associated with previously identified haploblocks showing strong differentiation between locally adapted populations. These extreme levels of standing genetic variation are likely associated with large effective population sizes and may help explain the remarkable adaptive divergence among populations of the Atlantic silverside.